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HALO  FORMATION  AND  HOLLOWING  IN 
RELATIVISTIC  ELECTRON  BEAMS 


I .  INTRODUCTION 

Many  of  the  analytical*  ^  and  computational^-^  models  used  to  treat 
relativistic  electron  beam  propagation  in  resistive  media,  such  as  air, 
assume  that  the  radial  profile  of  the  beam  current  density  J^(r)  has  a 
specified  shape.  The  Bennett  profile,  which  has  the  form, 


(1) 


is  often  assumed;  it  is  the  correct  equilibrium  for  a  beam  with  no  plasma 

current  and  an  isothermal  Maxwellian  distribution  of  perpendicular 
8  9  1 

velocities.  ’  An  envelope  equation  is  often  used  to  describe  the 

evolution  of  the  scale  radius  a^  as  the  beam  propagates. 

If  the  beam  induces  a  plasma  current  Jp(r),  the  Bennett  profile 

defines  the  equilibrium  only  if  Jp(r)  has  the  same  profile  as  the  beam. 

As  revealed  in  particle  simulation  code  studies,*0-*8  compact  Bennett-like 

beam  current  density  profiles  generally  occur  when  the  plasma  current  is 

small  or  when  its  profile  approximates  the  local  J^fr).  Such  codes  contain 

a  model  for  generating  the  conductivity  a(r)  as  the  beam  ionizes  the  air  or 

other  gas,  and  since  J  =  ctE  and  the  axial  electric  field  E  usually 

p  z  z 

varies  slowly  with  r,  the  plasma  current  profile  tends  to  follow  the 

4  5 

conductivity  profile.  ’  If  the  local  conductivity  is  either  much  narrower 

or  broader  than  the  beam,  and  if  the  plasma  current  Ip  is  a  significant 

12 

fraction  of  1^,  large  departures  from  the  Bennett  profile  may  occur. 

One  widely  studied  example  of  this  occurs  when  E is  large  enough  to 

cause  strong  avalanche  ionization  near  the  beam  axis. *°’ *^^0  Joyce  and 

Lampe*0  showed  that  this  process  results  in  violent  hollowing,  which  was 

first  observed  in  particle  simulation  codes,  and  later  confirmed 

20 

experimentally  by  Ekdahl,  et  al.  These  experiments  also  confirmed 
Manuscript  approved  December  7.  1987 
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theoretical  predictions  that  the  instability  could  be  eliminated  by 
limiting  the  rate  of  rise  of  the  beam  current  at  the  front  of  the  pulse. 
Although  the  hollowing  process  resembles  a  convective  instability,  it  is 
more  correct  to  state  that  the  beam's  original  Bennett-like  profile  ceases 
to  be  an  equilibrium,  and  the  beam  is  seeking  a  new  hollowed  equilibrium. 

A  different  process  occurs  when  the  conductivity  profile  is  broader 
than  the  beam.  Simulations  of  such  beams  often  evolve  to  a  "core  and  halo" 
equilibrium  in  which  part  of  the  beam  is  confined  to  a  small  region  near 

12 

the  axis  while  most  of  the  current  flows  in  a  halo  of  much  larger  radius. 

Even  if  the  plasma  current  fraction  is  small,  this  process  is  observed  to 

12 

occur  if  the  transverse  beam  temperature  is  too  large  at  injection.  Low 

21 

pressure  FX-100  propagation  experiments  by  Ekdahl  may  have  been  the  first 

to  exhibit  the  core -halo  phenomenon.  Cerenkov  data  taken  during 

propagation  experiments  with  the  RADLAC  accelerator  sometimes  show  what 

appears  to  be  a  halo  after  several  meters  of  propagation.  Perhaps  the 

strongest  evidence  for  core-halo  currents  comes  from  recent  x-ray  probe 

23 

measurements  made  during  low  pressure  propagation  experiments  on  ATA. 

This  paper  describes  analytical  and  simulation  studies  of  hollowing 

and  halo  formation  processes.  Analytical  models  are  developed  for- 

estimating  the  equilibrium  J^(r)  for  an  isothermal  beam,  in  the  presence  of 

a  plasma  current  with  a  specified  Bennett  profile  with  radius  ap.  A  purely 

analytical  solution  is  obtained  by  iterating  about  the  vector  potential 

generated  by  a  Bennett  beam  with  radius  a^Q.  On-axis  hollowing  is 

predicted  when  1^  +  1^  «  and  a^Q  is  somewhat  larger  than  a^.  Large 

radial  wings  ("halos")  are  generated  when  a^Q  is  small  compared  with  a^. 

In  extreme  cases  with  a,  <<  a  and  I,  +  I  <<  I.  ,  the  beam  current  profile 

bo  p  b  p  b 

J^(r)  develops  an  off-axis  minimum.  We  also  show  results  from  a  more  exact 

model  which  numerically  integrates  Ampere's  equation  as  described  by  Sharp, 
3 

et  al.;  it  shows  qualitatively  similar  behavior. 


2 


Later  portions  of  this  paper  describe  simulations  of  these  phenomena 

using  the  axisymmetric  code  SIMM0^  and  the  3-D  relativistic  simulation 
11-13 

code  SARLAC.  Since  the  on-axis  hollowing  process  has  been  widely 

studied,  most  of  the  results  are  for  cases  in  which  the  beam  evolves  into  a 
core  and  halo  equilibrium.  Halo  formation  is  most  common  in  high  current 
beams  (>  20  kA)  because  they  produce  a  large  return  current  fraction  and 
also  tend  to  have  broad  conductivity  profiles  because  of  recombination 
effects.^  Halo  formation  is  further  enhanced  if  the  beam  temperature  or 
emittance  at  injection  is  too  large.  Mismatched  beams  of  this  sort  can 
produce  substantial  halos  even  in  cases  with  modest  return  current 
fractions. 

The  last  section  discusses  the  effect  of  halo  equilibria  on  the 

5  5_7 

resistive  hose  instability.  Analytical  models  and  simulations  which 

assume  a  Bennett  profile  for  the  beam  predict  a  substantial  reduction  in 

hose  instability  growth  as  the  beam  current  increases.  Simulations  with 

the  SARLAC  code  generally  show  hose  growth  decreasing  as  1^  increases,  but 

much  more  weakly.  The  SARLAC  simulations  suggest  that  the  halo  is  a 

12 

destabilizing  factor,  since  it  contains  particles  whose  individual 
betatron  frequencies  can  resonate  with  low  frequency  hose  modes  produced  in 
the  expanded  head  of  the  beam.  Simulations  with  the  VIPER  code,^  which 
assumes  a  Bennett  profile  for  the  beam,  generally  agree  well  with  SARLAC  in 
cases  in  which  the  beam  remains  Bennett-like.  However,  VIPER  predicts  much 
less  instability  growth  than  does  SARLAC  in  high  current  or  poorly  matched 
cases  which  evolve  to  a  core  and  halo  equilibrium.  Since  halo  formation 
has  a  strong  effect  on  hose  stability,  it  deserves  more  experimental  study 
in  the  future. 
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II.  ANALYTIC  EQUILIBRIUM  MODELS 


2  3 

A.  Assumptions  and  formalism:  Lee,  and  Sharp,  Lampe  and  Uhm  ,  have 

discussed  isothermal  equilibria  for  ultrarelativistic  electron  beams 

propagating  through  dense  gases.  Conventionally,  the  beam  is  assumed  to  be 

axisymmetric,  and  large  aspect  ratio  approximations  are  made,  i.e., 

v^  <<  v^  =  c  for  all  particles  (the  "paraxial"  approximation), 

a^  3/3z  «  1,  and  (a^/c)  3/ 9 1  <<  1.  It  is  convenient  to  use  r,  z,  and 

C  =  ct-z  as  independent  coordinates  instead  of  r,  z,  and  t;  the  quantity  C 

is  the  distance  behind  the  beam  head,  and  is  a  constant  of  the  motion  if 

v^  =  c.  For  beam  injection  into  neutral  gas,  the  gas  conductivity  o(r,C,z) 

is  due  to  ionization  by  the  beam  itself.  It  is  assumed  that  a  is  large 

enough  to  short  out  radial  electric  fields  within  the  beam.  This  is 

generally  well  satisfied,  except  at  the  very  front  of  the  beam  and  at  large 
2  3 

radii.  ’  The  fields  Ez  and  are  then  derived  from  A,  the  z-component  of 
vector  potential: 

Ez  =  -  3A/3C, 

BQ  =  -  3A/3r, 

and  A  is  determined  by  Ampere’s  Law, 


1  3_ 
r  3r 


3A 

3r 


while  J  is  determined  by  Ohm's  Law, 
P 


J  =  -  a  3A/3C 
P 


(2a) 


(2b) 


We  use  the  conventions  that  J^r)  and  the  total  beam  current  are 

positive,  while  J ^ ( r )  and  the  total  plasma  current  1^  are  negative.  We  also 

define  a  (positive  constant)  return  current  fraction  f  =  -  I^/I^,  a  net 

current  I  s  I.  +  I  ,  and  a  net  current  fraction  f  s  I  /I.  s  1-f. 
n  b  p  n  n  b 

Ordinarily,  0  <  f  <1. 

J  '  n  - 
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In  our  analytic  modeling,  we  shall  assume  that  J^(r,0  is  C- 

independent  and  also  use  the  convenient  fiction  that  Jp(r,C)  is  independent 

of  £.  In  fact,  we  shall  assume  a  particular  form,  Eq.  (6)  below,  for 

J  (r).  Strictly  speaking,  Eqs.  (2)  specify  J  and  imply  that  J  must  vary 
P  P  P 

with  C,  but  this  variation  is  often  quite  slow  in  the  main  body  of  the 

3  5 

beam,  and  thus  it  is  quite  reasonable  to  assume  is  C-independent .  ’  In 
the  particle  simulations  discussed  in  Sec.  Ill,  J^(r,C),  o(r,0  and  Jp(r,C) 
are  all  determined  self-consistently,  and  we  do  not  assume  J^(C)  or  Jp(C) 
to  be  ^-independent . 

In  our  analytic  model,  the  beam  is  also  assumed  to  have  a  Maxwellian 
distribution  of  transverse  velocities  with  a  uniform  temperature  T  given 
by2'3 


-  1  -21  2  AQ  2 
T  *  2  Y  me  V1  £  2  Y  V  A9t 


(3) 


Here  v^  =  cA9t  is  the  transverse  thermal  velocity  of  the  beam.  This 

isothermal  Maxwellian  assumption  is  reasonable  for  a  beam  which  has  passed 

through  a  thick  scattering  foil  at  injection,  or  which  has  been  scattered 

9 

while  propagating  through  uniform  gas.  (In  the  simulations  of  Sec.  Ill, 
the  isothermal  Maxwellian  distribution  of  v^  is  taken  only  as  an  initial 
condition,  and  the  beam  is  free  to  evolve  to  a  different  distribution  if  it 
wishes.)  The  radial  force  on  a  beam  electron  is  -  e  3A/3r;  thus  in 
Maxwell-Bol tzmann  equilibrium  the  beam  current  density  Jb(r)  must  satisfy 


Vr> 

where  J. 

bo 

Jh(r)  has 
potential 


=  Jbo  exp  [ eA( r )/T J  (4) 

=  Jb<r=0).  Equations  (2)  and  (4)  are  satisfied  if  Jp  =  0  and 
the  Bennett  profile  defined  by  (1).  In  this  case,  the  vector 
is  given  by 


5 


(5) 


A(r)  =  -  -f  In 


lb  ' 


B.  Iterative  Solution  for  a  Non-Bennett  Beam:  An  iterative  solution  for 
J^(r)  can  be  obtained  for  an  isothermal  beam  with  a  Bennett  plasma  current 
profile  given  by 


(6) 


The  plasma  current  profile  is  assumed  to  be  fixed,  and  an  iteration  is 
performed  which  begins  with  a  zeroth  iterate  J^°^(r)  in  the  form  of  a 
Bennett  beam  profile  (1)  whose  radius  a^Q  can  be  larger  or  smaller  than  a^. 
This  is  used  in  Eq.  (2a)  to  calculate  the  zeroth  iterate  to  the  vector 
potential , 


A<»>(r)  -  -  ^  In  (l  •  ^  1"  ('  *  '  <7> 

This  is  substituted  into  (4)  to  obtain  the  next  iterate  to  the  equilibrium 
beam  current  density, 

-2T  / T  2fT„/T 

■  Jbot1  *  rSo2)  t1  *  r2/ap2)  •  <8> 


Here  T„  s  el, /2c  is  the  Bennett  temperature. 

It  is  important  to  note  that  a^  and  T  are  not  independent  quantities: 
a^  must  be  chosen  so  that  J^0  satisfies  the  Bennett  pinch  condition,  which 
is  a  necessary  but  not  sufficient  condition  for  equilibrium.  (The  Bennett 
pinch  condition  guarantees  that,  in  a  radially  averaged  sense,  the  magnetic 
pinch  balances  the  beam  pressure,  but  it  does  not  guarantee  that  pressure 
balance  holds  at  each  radial  location.)  This  condition  is 


6 


(9a) 


T  =  el  /2c, 
e  ’ 


where  the  "effective  current"*’^’'1  I  is  defined  by 


Ie  s  2Ib  dr  2 nr  Jbv°;(r)  Inv°;(r), 


and  for  the  particular  case  where  Jb(r)  and  J p ( r )  are  given  by  Eqs.  (1) 
and  (6), 


I  =  21, 

e  b 


dr  2 nr 


-  2  h  2/  2)2  1,  2,  2  T  .  2. 

nabo  l1  +  r  /abo  J  U  +  r  /abo  1  +  r  /a 


2/  2)2 
r  /abo  J 


2,  2 
i  1  /a, 
b _ bo 

2,  2 
+  r  /a, 
bo 


r  2.  2 

I  r  /a 


2.  2 

/  a 

P 


=  *b  1 


2,2 
2fa,  /a 
bo  p 


ho2'^2)-1 


f,  ln(ab°;/aPZ) 

)->l  kV)-1 


Here  1^  7(r)  is  the  net  current  enclosed  within  radius  r.  is 

essentially  a  measure  of  the  pinch-ef fectiveness  of  the  net  current 

contained  within  the  beam,  whereas  I  s  I,  +  1  is  the  total  net  current 

- - —  -  n  b  p 

out  to  r  ■  “  (or  to  the  walls).  Thus,  T  >1  if  a,  <  a  (so  that  much  of 

e  n  bo  p 

the  plasma  return  current  is  outside  the  beam).  1^  *  Jb  in  the  limit 

a,  /a  -»  0.  Conversely,  I  <  I  if  a,  >  a  ,  and  I  -»  I,  +  21  when  a,  /a 
bo  p  e  n  bo  p  e  b  p  bo  p 

->  m.  This  is  because  current  elements  that  are  entirely  within  the  beam 

are  twice  as  effective  in  determining  the  average  pinch  force.  Thus  one 

may  regard  Eqs.  (9)  as  a  relation  determining  abo  as  a  function  of  T  and  f, 

or  alternatively  defining  T  as  a  function  of  a,  /a  and  f.  We  note  that 

bo  p 

increasing  T  corresponds  to  reducing  abo  (in  order  to  minimize  the  plasma 
return  current  contained  within  the  beam).  This  is  at  odds  with  intuition, 
since  one  usually  thinks  of  increased  beam  transverse  temperature  as 


leading,  in  a  dynamic  sense,  to  expansion.  We  note  that  no  equilibrium  is; 
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possible  {(9)  has  no  solution)  if  T  >  Tg  or  T  <  Tg(l-2f).  In  the  former 
case,  the  beam  is  too  hot  to  be  pinched  even  by  the  unshielded  beam  pinch 
force,  whereas  in  the  latter  case  the  beam  is  too  cool  to  support  itself 
even  against  the  weakened  pinch  force  corresponding  to  return  current  flow 


entirely  in  a  narrow  region  on-axis.  Figure  1  plots  T  versus  a^/a^  for  a 

series  of  cases  with  f  =  0.15  and  f  =  0.8. 

n  n 

Equations  (8)  and  (9)  give  a  great  deal  of  detailed  information  about 


the  beam  profiles  J^(r).  This  information  is  summarized  in  Table  I  and 
will  be  discussed  in  this  subsection  and  the  next. 


In  cases  with  small  return  currents,  the  equilibrium  is  expected  to 

remain  close  to  the  Bennett  profile,  which  is  the  exact  solution  if  I  =  0. 

P 

Figure  2  plots  profiles  of  ^ (r )/Jg^ ^ ( r=0)  from  (8)  using  ab0/ap  =  0.4, 

0.7,  1.0,  and  1.15,  and  T  calculated  from  (9).  The  net  current  fraction  is 

n>  -4 

£n  =  0.8.  In  all  cases,  Jg  (r)  is  close  to  the  r  scaling  expected  for 
a  Bennett  beam  at  large  radius;  no  hollowing  or  halo  formation  is  observed. 

Next  we  consider  cases  in  which  the  return  current  is  substantial, 
beginning  with  the  regime  a^0/ap  <  1*  for  which  1-f  <  T/Tg  <  1.  Using 
Eqs.  (9),  this  condition  can  also  be  written  in  the  equivalent  form 
f  <  f  ,  which  indicates  directly  that  in  some  average  sense  the  beam  is 
narrower  than  the  Jp(r)  profile.  Here,  Eq.  (8)  shows  that  in  the  central 
part  of  the  beam,  J^^(r)  scales  as  (1  +  r^/a^^)  P  wi  th  p  >  2,  and  thus 
has  a  steeper  fall-off  than  a  Bennett.  However,  J^^(r)  scales  as  r  ^  , 
with  p'  <  4,  for  large  r.  Thus,  there  is  a  tightly  pinched  central  core 
sutrounded  by  an  extended  lower  density  halo.  It  is  interesting  to  note 
that  as  T  increases  above  (l-f)Tg,  the  core  radius  decreases  (as  predicted 
above  in  the  discussion  of  ) ,  but  progressively  more  of  the  current  is 
carried  in  the  halo,  whose  radius  increases .  Thus,  the  half  radius 


(defined  as  the  radius  within  which  half  the  beam  current  flows)  typically 
increases  as  a  function  of  T. 
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Figure  3  shows  a  series  of  calculations  for  a  case  with  a  strong 
plasma  return  current,  fn  =  0.15.  The  beam  equilibrium  for  a^/a^  =  0.4 
exhibits  core  and  halo  behavior:  J^^(r)  drops  by  three  orders  of 
magnitude  between  0  <  r/a^  <  3,  but  only  by  a  factor  of  3  as  r/a^  increases 
from  3  to  10. 

When  a^/a^  =  1,  it  follows  that  T  =  Tg(l-f)  and  Eq.  (8)  shows  that 
the  Bennett  profile  is  recovered.  This  is  shown  as  the  second  curve  in 
Fig.  3. 

For  a^/a^  >1,  Eq.  (8)  gives  two  subcases.  If  J^q(0)  +  J^fO)  <  0,  the 
profile  J^^(r)  is  hollowed  on-axis,  i.e.,  has  a  minimum  on-axis  and  a 
single  maximum  somewhere  off-axis.  This  can  be  seen  by  differentiating  Eq. 
(8);  extrema  are  found  at  r  =  0  and  at 


y0)  +  jb(o)(0) 


i.  + 1 

b  p 


2  2 
a,  a 
bo  p 


One  and  only  one  extremum  for  r  >  0  exists  if  and  only  if  -J^(0)  > 
Jb(0)((»,  which  can  occur  only  if  the  beam  is  wider  than  the  return  current 
profile.  This  situation  is  shown  in  the  top  curve  of  Fig.  3,  where 


a.  /a  =1.15  and  f  =  0.15.  If,  on  the  other  hand,  a,  /a  >1  but 
bo  p  n  bo  p 

Jb(o)0»  +■  Jp(0)  >  0,  then  J^^(r)  is  flattened  near  r  =  0,  but  is  not 
hollowed. 


The  equilibria  shown  in  Fig.  3  resemble  the  "generalized  Bennett" 
equilibria  obtained  by  Benford,  Book,  and  Sudan. ^  However,  their 
generalized  Bennett  J^fr)  is  zero  on  axis  for  all  cases  which  exhibit  on 
axis  hollowing,  and  their  cases  which  resemble  core/halo  profiles  have 
infinite  current  density  at  r  =  0. 
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Finally,  we  consider  the  occurrence  of  off-axis  hollowing,  which  is 
defined  as  having  a  ( r )  profile  that  has  a  maximum  on-axis,  a  minimum 

somewhere  off-axis,  and  then  a  second  maximum  further  out  radially.  Since 


Eq.  (8)  gives,  at  most,  one  off-axis  extremum,  off-axis  hollowing  is  not 
predicted  by  the  first  iteration  J^^(r).  However,  a  second  iteration 


does  predict  off-axis  hollowing  under  conditions  where  is  small  and 


f  is  small.  We  can  see  this  by  going  part  way  towards  a  next  iteration. 


In  equilibrium,  the  beam  satisfies  the  pressure  balance  condition 


1  t  3P 

c  bz  0  3r 


_  L-  61) 

3r  lec  J 


Hollowing  occurs  when  there  is  a  region  where  3J^/3r  >  0.  According  to 
(11),  this  requires  that  Bg(r)  change  sign,  which  is  equivalent  to 
requiring  that  the  enclosed  net  current  Ip(t)  change  sign.  To  first  order, 
the  net  current  density  Jn^(r)  -  Jb^\r)  +  Jp(r)  is  given  by  (6)  and 
(8),  and  the  enclosed  net  current  is 

In<1)(r)  =  Ib(1)<r)  +  I  <r) 


=  2 rt  dr  r  1  + 


212fTB/T 


~,-2TD/T  _  2.  2 

2]  B  I  r  /a 


,  2.2 

1  +  r  /a 

P 


The  integral  in  (12)  can  be  evaluated  numerically. 

Figure  4  plots  the  enclosed  currents  I^^(r),  Ip(r)  and  In^\r)  for 

a  case  with  a,  /a  =  0.8  and  f  =  0.15.  The  enclosed  net  current  remains 
bo  p  n 

positive  everywhere,  so  this  case  is  expected  to  show  core  and  halo 
behavior.  In  Fig.  5,  abQ/ap  has  been  reduced  to  0.4,  and  there  is  a 
significant  region  where  I  ^\r)  is  negative.  This  example  would  thus 
show  off-axis  hollowing  in  the  next  iteration. 


10 


C.  Numerical  Integration  of  Ampere's  Equation:  The  iterative  procedure 
offers  considerable  insight  into  the  hollowing  and  halo  formation 
processes,  but  one  or  two  iterations  does  not  necessarily  yield  a 
quantitatively  accurate  calculation  of  Jb(r)  when  ( r )  is  very  different 
from  the  Bennett  profile  which  is  the  starting  point  for  the  iteration.  In 
the  present  section  we  continue  to  assume  that  Jp(r)  has  the  Bennett  shape 
(6),  but  we  calculate  Jb(r)  directly  by  numerically  integrating  Ampere's 
law  (2a)  and  the  Boltzmann  equilibrium  condition  (4)  together.  This  method 

3 

was  employed  by  Sharp,  et  al.,  to  construct  equilibria  for  hose 
instability  calculations.  Values  for  a^,  Jp(0),  J ^ ( 0 )  and  T  are  chosen, 
and  (2a)  is  integrated  numerically  from  r  =  0  to  f  =  %  using  (4)  to 
generate  Jb(r).  Ue  then  calculate  by  integrating  Jb(r): 


03 


-  2n 


dr 


r  Jb(r) 


(13) 


Jo 

and  f  s  -  I  /I.  or  f  s  1-f  are  thus  specified.  The  procedure  can  be 
p  b  n 

repeated,  adjusting  the  initial  value  of  Jb(0)/Jp(0),  to  converge  to  a 
desired  value  of  f,  or  alternatively  to  a  desired  ratio  of  a^  to  some 
characteristic  beam  radius. 

Figures  6-8  show  examples  in  which  fn  <  T/Tg,  and  thus  the  core  of  the 

beam  has  a  narrower  profile  than  Jp(r)  and  beam  profiles  of  core/halo  type 

are  expected.  In  Fig.  6,  f  =  0.47  and  T/Tg  -  0.6;  because  the  return 

current  is  not  too  large  Jb(r)  is  close  to  the  Bennett  profile,  with  a  weak 

halo.  In  Fig.  7,  T/Tg  =  0.2  and  there  is  a  stronger  return  current, 

f  =  0.18.  This  beam  more  clearly  exhibits  core  and  halo  behavior,  as 
n 

expected  from  the  simple  iterative  theory  of  Sec.  IIB.  However,  the 
present  solution  (in  contrast  to  the  more  approximate  result  shown  in 
Fig.  2)  shows  that  the  halo  does  not  extend  indefinitely,  but  rather  Jb(r) 


falls  off  rapidly  for  r  >  2  a^.  This  type  of  behavior  is  seen  in  the 
simulations  reported  in  Sec.  III.  In  Fig.  8,  T/T„  =  0.1  and  the  return 

D 

current  is  still  stronger,  f  =  0.083.  This  example  exhibits  off-axis 
hollowing.  Figure  9  shows  a  case  in  which  there  is  strong  return  current, 
fn  =  0.08,  and  T/Tfi  =  0.05  <  fn>  Thus  Jp(r)  has  a  narrower  profile  than 
the  beam,  and  in  fact  -Jp(0)  >  J^(0) -  The  beam  exhibits  on-axis  hollowing, 
as  expected. 
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III.  PARTICLE  SIMULATION  RESULTS 


A.  Model  Descriptions:  Particle  simulations  have  been  widely  used  to 

9-17 

study  relativistic  electron  beam  propagation.  Most  of  the  results 

presented  here  were  obtained  with  SARLAC,  an  ultrarelativistic  3-D  particle 
11-13 

simulation.  It  differs  from  conventional  plasma  simulations  in  that 

v^  «  vz  =  c,  so  particles  remain  at  the  same  distance  C  from  the  beam 

head,  and  field  and  conductivity  equations  are  integrated  in  C*  The  field 

5  25 

solver  is  based  on  the  Lee  field  equations  ’  for  a  relativistic,  paraxial 

26 

beam.  The  conductivity  equations  employ  a  variable  temperature  model  and 

include  beam  impact  ionization,  avalanche,  and  recombination.  The  beam  is 

injected  with  a  Bennett  profile  with  radius  aQ  and  a  constant  temperature 

11  13 

T  .  The  SARLAC  code  is  described  in  more  detail  elsewhere.  ’ 
o 

Results  on  on-axis  hollowing  were  obtained  with  SIMM0,  an  axisymmetric 
particle  simulation. ^  Except  for  the  axisymmetric  assumption,  the  code  is 
similar  to  SARLAC. 


B.  Overview  of  SARLAC  Simulations:  Eight  simulations  were  performed  with 
the  SARLAC  code  using  various  values  of  beam  current  and  beam  transverse 
temperature.  These  simulations  are  described  in  Table  II.  The  high 
current  (>  100  kA)  cases  produced  large  return  currents  and  broad 
conductivity  profiles  (due  to  recombination  effects  and  beam  tapering).  In 
agreement  with  the  theory  presented  in  Sec.  II,  these  cases  show  core  and 
halo  behavior  when  the  initial  transverse  temperature  Tq  is  large,  i.e., 


f 

m 


s  T  /Tn  > 
o  B 


f  . 
n 


(14) 


In  simulation  parlance,  f  defined  in  (14)  is  called  the  "matching  current 

fraction"  --  the  beam  as  injected  would  be  in  a  Bennet t- like  equilibrium  if 

the  net  current  fraction  f  turned  out  to  be  equal  to  f  .  If  the  "matching 

n  ^  m  6 
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current"  is  larger  than  the  self-consistently  determined  net  current,  that 
means  that  the  beam  is  "overpressured"  -  i.e.,  too  hot  -  at  injection. 

Halos  were  produced  in  a  10  kA  case  when  £  was  very  large,  even  though  the 

!T1 

return  current  fraction  was  modest.  In  one  extreme  case  with  1^  =  500  kA, 
off-axis  hollowing  was  observed. 

To  quantify  the  degree  of  halo  formation,  the  radius  a^/2  enclosing 

half  of  the  beam  current  was  compared  with  a  ,  a  Bennett  radius  fit  to 

core 

the  beam  core.  For  a  Bennett-like  beam,  a,  =  a _ ,  while  a.  >2  aD 

1/z  core  l/z  -  B 

for  a  beam  with  a  strong  halo.  The  radius  measurements  were  taken  at 

C  =  80  a  and  z  =  240  a  where  a  was  the  Bennett  radius  of  the  injected 
ooo 

beam. 


C.  Bennett-like  Examples:  Case  1  is  a  10  kA  beam  with  f  =  0.5,  a 

temperature  close  to  that  required  for  a  well-matched  injection.  Since 

Ip/I^  =0.5  for  this  case,  the  beam  should  remain  Bennett-like.  Figure  10 

plots  J . ( r )  and  J  (r)  s  J.(r)  +  J  (r)  at  C  =  80  cm  and  z  =  240  cm.  The 
r  b  n  b  P 

profiles  are  similar  in  shape  and  there  is  no  evidence  for  a  halo.  Figure 

11  compares  J^(r)  with  a  Bennett  fit;  the  fit  is  excellent  even  out  to 
large  radii  where  ( r )  has  dropped  by  several  orders  of  magnitude.  Figure 

12  plots  the  pinch  force  per  beam  electron  F^fr)  for  this  case.  It  shows 
the  expected  smooth  behavior  for  a  Bennett-like  beam. 

Figure  13  plots  the  x  and  y  positions  of  a  subset  of  individual 
simulation  particles  vs  C  after  the  beam  has  propagated  to  z  =  320  cm.  The 
beam  flares  at  the  head  because  it  takes  time  for  the  net  current  and  pinch 
force  to  build  up.^  The  beam  is  clearly  well-pinched,  and  hose  instability 
distortions  are  barely  visible. 
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A  100  kA  beam  typically  produces  a  broad  conductivity  profile  and  a 
small  net  current  (In  ~  0.1  -  0.2  1^),  conditions  which  typically  lead  to  a 
core/halo  equilibrium.  SARLAC  simulations  of  100  kA  beams  reached  Bennett- 
like  equilibria  only  when  the  head  of  the  beam  was  expanded  and  cooled 
prior  to  injection  into  the  gas,  in  order  to  improve  the  match  throughout 
the  beam.  In  cases  6  and  7,  the  beam  was  injected  with  a  tapered  radius 

a( C»  z=0)  =  (tr/t)1/2ao  ,  (15) 

but  the  emittance  at  injection,  which  is  proportional  to 
1/2 

a(C,0)(f  (t,0))  ,  was  constant  as  a  function  of  C.  For  f  =  0.2  (case 

m  m 

7),  the  beam  remains  Bennett-like,  as  shown  in  Fig.  14.  For  f  =  0.4  (case 
6),  modest  departures  (core/halo)  from  the  initial  Bennett  profile  occur. 

D.  Core  and  Halo  Examples:  Case  4  is  a  100  kA  beam  with  f^  =  0.4,  a 

relatively  high  temperature  for  a  beam  with  such  a  small  net  current 

fraction.  Figure  15  plots  J^r)  and  Jn(r),  showing  the  presence  of  a  halo 

between  r  =  5  cm  and  r  =  15  cm.  The  beam  is  well  represented  by  a  Bennett 

out  to  r  =  5  cm.  Since  a^/ ^  >  3acoi_g  (Table  II),  this  may  be  considered  a 

strong  halo  case.  The  core,  whose  radius  a  =  1.1  a  =3.3  cm,  carries 

core  o 

only  about  15%  of  the  total  current.  Most  of  the  beam  is  contained  in  the 
halo.  The  net  current  density  Jn(r)  is  actually  negative  just  inside  the 
halo.  This  leads  to  a  local  off-axis  minimum  in  | ( r )  | ,  as  seen  in 
Fig.  16.  This  tends  to  push  particles  from  the  core  to  the  halo.  The  halo 
is  clearly  seen  in  the  particle  plots  shown  in  Fig.  17,  along  with  the  high 
density  core.  Consistent  with  Fig.  15,  the  halo  has  a  relatively  well- 
defined  edge.  Codes  such  as  VIPER,  which  impose  a  Bennett  profile  on  the 
beam  and  use  an  envelope  equation  to  estimate  the  beam  radius,  predict  that 
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the  beam  radius  should  taper  back  from  the  beam  head  due  to  the  gradual 
build-up  of  the  effective  current  defined  in  (9).  The  behavior  seen  in 
Fig.  17  is  much  different.  The  implications  of  this  for  hose  instability 
are  discussed  in  Sec.  IV. 

Strong  halo  formation  of  a  similar  nature  is  seen  in  Case  3 

(I,  =  50  kA,  f  =  0.4).  In  Case  5  (I,  =  100  kA)  f  was  reduced  from  0.4  to 
b  m  b  m 

0.2  in  order  to  produce  a  better  match  at  injection;  only  mild  halo 

formation  was  observed  in  this  case. 

We  have  also  observed  halo  formation  in  lower  current  beams  when  the 

initial  temperature  or  emittance  is  poorly  matched.  Case  2  is  such  an 

example;  it  is  identical  to  the  10  kA  Case  1  except  that  f  =  1.5.  As  seen 

in  Fig.  18,  a  halo  forms  at  r  =  2  cm,  and  the  Bennett  core  carries  only 

about  25%  of  the  total  current.  From  Table  II,  a,  ,,/a  =  2.5  for  this 

1/2  core 

case.  We  believe  that  this  case  is  qualitatively  different  from  the  high- 
current  cases  discussed  in  this  paper;  there  is  not  enough  plasma  return 
current  to  account  for  the  halo  through  the  mechanism  discussed  in  Sec.  II. 
In  this  case,  the  evolution  to  a  core/halo  equilibrium  appears  to  be  due  to 
a  strong  departure  from  Maxwell-Bol tzmann  equilibrium.  This  beam  is 
injected  with  a  transverse  pressure  that  far  exceeds  the  pinch  force,  and 
the  resulting  non-adiabatic  behavior  (large-amplitude  bounces)  heats  some 
electrons  far  more  than  others  and  introduces  structure  in  phase  space.  We 
believe  that  this  type  of  behavior  may  be  a  significant  feature  of  some 
beam  propagation  experiments,  and  intend  to  explore  it  further  in  future 
work. 


E.  Off-Axis  Hollowing  Cases;  Figure  8  indicates  that  off-axis  hollowing 

can  occur  if  f  s  I  /I,  is  sufficiently  small.  Small  f  is  associated  with 
n  n  b  n 

very  high  beam  current.  Case  8  is  similar  to  Case  2  but  with  I,  =  500  kA 

b 
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and  fm  =  0.2.  A  plot  of  J^r)  at  z  =  1200  cm  and  £  =  480  cm  (Fig.  19) 

reveals  an  off-axis  minimum  at  r  =  40  cm.  Again  the  core  contains  only  a 

small  fraction  of  the  current.  An  off-axis  minimum  is  not  observable  in 

plots  of  the  x  and  y  positions  of  simulation  particles,  so  Fig.  20  also 

plots  the  radial  position  r^  of  simulation  particles  at  z  =  200  aQ  =  1200 

cm.  The  core  has  a  relatively  sharp  edge,  followed  by  a  gap  in  particles 

corresponding  to  the  off-axis  minimum. 

The  conditions  to  produce  off-axis  hollowing  by  a  beam  injected  into 

an  initially  neutral  gas  are  well  beyond  present  technology.  However,  our 

axisymmetric  simulations  using  the  SIMM0  particle  code  indicate  that  off- 

axis  hollowing  can  occur  when  a  lower  current  beam  is  injected  into  a 

highly  conducting  channel,  which  carries  a  current  1^  ~  0.1  1^  and  is 

somewhat  broader  than  the  beam.  Such  conditions  can  occur  when  a  beam  is 

27 

injected  into  a  laser-guided  electric  discharge.  Keeley  has  observed 
similar  behavior  using  a  fluid  beam  simulation  code. 

F.  On-Axis  Hollowing:  As  seen  in  previous  theoretical  and  experimental 
studies,  on-axis  hollowing  usually  causes  violent  disruption  of  the 
propagating  beam.^’^  The  disruptive  instability  is  driven  by 
enhancement  of  avalanche  ionization  near  the  beam  axis.  An  example  of  on- 
axis  hollowing  which  remains  quiescent  and  near  equilibrium  is  shown 
in  Fig.  21.  The  simulation  was  performed  with  the  SIMM0  code  for 
a  10  kA  beam  propagating  in  a  reduced  density  channel  with  density 

p(r=0)  =  0.1  atm.  The  channel  radius  a  is  twice  the  beam  radius,  and 

c 

p(r  >>  ac)  approaches  1  atm.  The  surface  plot  of  J^(r,C)»  at  z  =  800  cm, 
exhibits  the  mild  on-axis  hollowing.  Figure  22  plots  the  radius 
al/2<C»  z).  Although  some  oscillations  are  present,  the  beam  propagates 
without  disruption. 
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IV.  HOSE  INSTABILITY  SIMULATIONS 


A.  Introduction:  The  resistive  hose  instability  can  seriously  disrupt 

2  3  5-7 

propagation  of  relativistic  electron  beams.  '  ’  However,  it  is 
sometimes  possible  to  avoid  severe  disruption  by  taking  advantage  of  the 
convective  nature  of  the  instability.  An  initial  excitation  convects 
backward  in  the  beam,  growing  as  in  convects.  Thus  the  hose  amplitude,  at 
a  given  point  C  in  the  beam,  initially  grows  and  then  reaches  a  maximum  and 
subsequently  decays  as  the  disturbance  convects  past;  this  is  known  as 
convective  saturation.  The  linear  theory  of  the  instability  has  been 
developed  in  Refs.  2,  3  and  5. 

SARLAC  was  primarily  intended  for  studying  this  instability.  It  is 
less  restrictive  than  older  phenomenological  codes  such  as  VIPER, ^ 

EMPULSE,^  and  PHLAP^  in  three  ways:  First,  SARLAC  can  follow  the 
instability  into  the  large-amplitude  regime,  whereas  the  older  codes 
linearize  in  hose  amplitude  (compared  to  beam  radius).  Second,  SARLAC 

follows  exact  particle  trajectories,  whereas  the  older  codes  use  the  spread 

2  3  5 

mass  or  multicomponent  beam  dynamics  models.  However,  careful  studies 

have  shown  that  these  models  are  quite  accurate,  at  least  for  small 

amplitude  hose.  Third,  the  older  codes  impose  a  Bennett  radial  profile  on 

the  beam,  but  SARLAC  allows  the  beam  profile  J ^ ( r )  to  change  its  shape 

self-consistently .  For  our  purposes  here,  this  is  the  most  important 

difference  between  SARLAC  and  the  older  codes.  In  this  section,  we  examine 

the  hose  instability  behavior  of  the  cases  described  in  the  previous 

section  and  compare  the  SARLAC  results  with  those  obtained  from  VIPER. 

B.  Bennett-like  Cases:  For  cases  in  which  the  beam  remains  Bennett-like, 
we  expect  VIPER  and  SARLAC  to  predict  comparable  hose  growth  in  the  linear 
regime.  Figure  23  plots  the  displacement  X(z)  of  the  beam  centroid  for 
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Case  1,  which  is  a  Bennett-like  10  kA  beam.  SARLAC  results  are  shown  for 

C  =  40,  80,  and  160  cm,  and  a  VIPER  result  for  C  =  160  cm.  (The  VIPER 

_2 

amplitudes  are  normalized  to  Xq  =  10  cm,  the  amplitude  of  the  initial 

SARLAC  perturbation).  As  expected,  the  two  codes  are  in  substantial 

agreement.  This  is  also  shown  in  line  1  of  Table  III,  which  compares  the 

SARLAC  and  VIPER  values  for  ln(X/X  )  evaluated  at  C  =  80  a  and 

o  o 

z  <  400  aQ.  Once  the  beam  propagates  a  short  distance,  the  radius 

al / 2 ( *n  t*le  re£i°n  strong  hose  growth  is  almost  constant.  This  is 

seen  in  Fig.  24. 

C.  Core  and  Halo  Examples:  In  all  cases  where  there  is  a  halo,  SARLAC 
predicts  more  hose  growth  than  VIPER.  In  Case  5,  only  mild  halo  formation 
is  seen  in  the  SARLAC  simulation,  and  SARLAC  shows  5.9  e-foldings  of 
growth,  a  mild  enhancement  compared  to  the  4.4  e- foldings  shown  by  VIPER. 
The  high  current  Case  4  exhibits  strong  halo  formation.  Figure  25  plots 
X(z)  for  this  case.  The  three  SARLAC  curves  at  C  =  120,  240,  and  360  cm 
show  substantial  instability  growth,  which  is  also  apparent  in  the  particle 
plots  shown  in  Fig.  26.  (This  plot  is  at  z  -  1200  cm  which  was  the  end  of 
the  simulation  run.  At  z  =  720  cm  (Fig.  17),  the  beam  does  not  yet  exhibit 
large  amplitude  hose. |  However,  the  normalized  amplitude  from  VIPER  at 

C  =  360  cm  shows  very  little  growth.  Indeed,  Table  III  shows  3.9  e-folds 
of  growth  for  SARLAC,  but  only  0.8  e-folds  for  VIPER.  Case  2,  a  10  kA 
beam  which  exhibits  strong  halo  formation,  also  shows  a  factoi  of  2  more 
instability  e-foldings  with  SARLAC  than  with  VIPER. 

D.  Discussion:  It  appears  to  be  generally  true  that  the  presence  of  a 
halo  exacerbates  hose  instability,  as  compared  to  the  same  case  treated 
with  a  model  that  artificially  imposes  a  Bennett  profile.  We  believe  this 
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is  true  because  J^r.C)  is  nearly  constant  in  both  r  and  C  within  the  halos 
seen  in  SARLAC  simulations,  and  thus,  there  is  much  less  variation  in 
betatron  frequency  among  beam  particles,  as  compared  to  beams  with  Bennett 
profiles.  It  has  been  known  for  many  years  ’  that  betatron  frequency 
spread  greatly  reduces  the  hose  growth  rate  and  leads  to  convective 
saturation.  Particularly  in  the  case  of  high-current  beams  (1^  »  10  kA), 
envelope  code  simulations  (which  impose  a  Bennett  profile)  typically  show  a 
gradual  tapering  of  the  half-radius  a^2  (0»  as  a  function  of  t,  as 
illustrated  in  the  VIPER  simulation  of  Case  4  in  Fig.  27.  This  tapering 
leads  to  betatron-frequency  variation,  and  is  a  stabilizing  feature.  The 
SARLAC  simulation  of  the  same  case  (Fig.  27)  shows  much  less  variation  in 
a^22  within  the  beam  body. 

It  should  be  noted  that  the  point  of  Table  III  is  to  compare  SARLAC  to 

VIPER  hose  growth  for  each  individual  case.  The  comparison  of  hose  growth 

among  different  cases  is  not  meaningful,  without  considering  many  factors 

that  have  not  been  discussed  here  and  that  are  different  for  the  various 

cases,  e.g.,  propagation  distances  z  must  be  scaled  to  the  betatron 

wavelength,  locations  C  (relative  to  the  beam  head)  must  be  scaled  to 

2 

"dipole  decay  lengths"  /2c,  and  the  variation  a^^^)  must  be  taken 

into  account. 
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V.  CONCLUSIONS 


Particle  simulations  of  propagating  relativistic  electron  beams  with 
large  return  current  fractions  often  exhibit  large  departures  from  the 
compact  Bennett  or  Gaussian  profiles  usually  assumed  in  theoretical 
calculations.  Beams  with  narrow  return  current  profiles  may  show  on-axis 
hollowing,  while  those  with  broad  return  current  profiles  may  form  a  core 
and  halo  or  an  off-axis  minimum  in  J^(r).  All  of  these  features  can  be 
reproduced  in  analytic  models  which  calculate  an  equilibrium  J,(r) 
compatible  with  Ampere's  equation. 

Beams  which  evolve  to  highly  non-Bennett  radial  profiles  may  be 
strongly  destabilized.  A  beam  with  a  narrow  conductivity  profile  and  an 
equilibrium  with  an  on-axis  minimum  is  likely  to  be  disrupted  by  the 
axisymmetric  hollowing  instability.  High  return  current  beams  with  broad 
conductivity  profiles  may  develop  a  halo  which  can  strongly  enhance  hose 
instability  disruption.  Since  most  of  the  beam  usually  resides  in  t he  halo 
in  these  cases,  the  presence  of  a  small  current  in  a  tight  core  may  be  of 
little  consequence.  The  important  role  of  the  J^(r)  profile  in  determining 
beam  stability  may  warrant  more  careful  experimental  examination  of  the 
phenomena . 
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Table  I.  A  summary  of  the  analytic  model  for  beam  equilibrium. 


f 

a.  /a 
bo  p 

t/tb 

Form  of  profile  J ^ ( r )  | 

0 

- 

1 

Bennett 

close  to  0 

any  value 

close  to  1 

Bennett-like 

>  1 

no  solution 

moderate 

<  1 

1-f  <  T/TB  <  1 

core/halo 

close  to  1 

<  1 

1-f  <  T/Tg  <  1 

off-axis  hollowed 

any  value 

1 

1-f 

Bennett 

s  2  .  2 
>  a  /a. 
p  bo 

>  1 

1 - 2 f  <  T/Tb  <  1-f 

hollowed  on  axis 

✓  2.  2 
<  a  /a. 
p  bo 

>  1 

1 - 2 f  <  T/TB  <  1-f 

flattened  on  axis 

<  ]-2f 

no  solution 
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Table  II.  This  table  summarizes  the  initial  conditions  and  the  type  of 

equilibrium  attained  in  eight  SARLAC  simulations.  In  all  cases,  the  beam 

electrons  initially  have  relativistic  factor  y  =  100,  the  gas  is  air  at 

standard  temperature  and  pressure,  the  simulation  region  is  terminated  by  a 

perfectly  conducting  wall  at  r  =  81  aQ,  and  the  beam  current  1^(0  rises  to 

its  full  value  over  a  rise  length  Cr  *  120  aQ.  Those  cases  designated 

"pencil  beam"  are  injected  with  uniform  emittance  and  radius  aQ,  while 

those  designated  "tapered"  are  injected  with  uniform  emittance,  but  with  a 

1/9  _1/9 

tapered  radius  a(C,)/ao  =  ,  and  T(C)  ~  [a(Q)  •  a1/2  is  the 

"half-radius"  of  the  beam,  while  a  is  a  Bennett  radius  fit  to  the 

core 

central  part  of  the  beam.  The  quantities  f  ,  a,  and  a  are  evaluated 

M  n’  1/2  core 

at  C  =  80  a  and  z  =  240  a  .  The  matching  current  fraction  f  is  defined 
o  o  e  m 

in  (14). 


Case 

Ib(kA) 

a  (cm) 

0 

f 

m 

£ 

n 

a,  ,,,/a 
1/2  o 

a  /a 

core  o 

Description 

1 

10 

1  (pencil) 

0.5 

0.53 

1.2 

i.i 

Bennett-like 

2 

10 

1  (pencil) 

1.5 

0.58 

3.5 

1.4 

Strong  halo 

1 

50 

3  (pencil) 

0.4 

0.21 

3.1 

1.5 

Strong  halo 

4 

100 

3  (pencil) 

0.4 

0.22 

3.5 

1.1 

Strong  halo 

5 

100 

3  (pencil) 

0.2 

0.16 

1.6 

0.9 

Moderate  halo 

6 

100 

3  (tapeLed) 

0.4 

0.17 

2.4 

1.7 

Moderate  halo 

7 

100 

3  (tapered) 

0.2 

0.13 

1.5 

1.8 

Bennet  t-like 

8 

500 

6  (pencil) 

0.2 

0.04 

3.2 

1.2 

Off-axis  hollowed 
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Table  III.  Hose  instability  growth,  for  four  of  the  cases  described  in  Table 


k 

t 

k 

f 

► 


2.  X  is  the  largest  hose  amplitude  seem  at  £  =  80  aQ  for  z  <  400  while 
Xq  is  the  initial  hose  amplitude. 


Case 

Ib(kA) 

Description 

SARLAC 

ln(X/Xo) 

VIPER 

ln(X/XQ) 

1 

10 

Bennett-lik.e 

5.4 

4.9 

2 

10 

Strong  halo 

0.6 

0.3 

4 

100 

Strong  halo 

3.9 

0.8 

5 

100 

Moderate  halo 

5.9 

4.4 
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0  0.2  0.4  0.6  0.8  1.0  1.2 

ab0/ap 

Fig.  1  Equilibrium  value  of  T/Tn  versus  the  nominal  beam  radius  ratio 

o 

a,  /a  .  The  solid  line  is  for  f  =  0.8,  and  the  dashed  line  is  for 
bo  p  n 

fn  =  0.15.  The  asterisks  correspond  to  the  cases  plotted  in  Figs. 

2  and  3. 
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Equilibrium  beam  current  density  ;(r)  calculated  from  (8)  for 

beam  with  net  current  fraction  f  =  0.8.  The  curves  are  for 

n 

nominal  beam  radius  a^  =  0.4,  0.7,  1.0,  and  1.15  (normalized  to 
the  plasma  current  radius  ap).  The  radial  profiles  are  all 


Bennett -like. 


I 


Fig.  4  Enclosed  currents  1^'  J(r),  Ip(r)  and  In'  ;(r)  versus  radius  i/a^ 

normalized  to  I^(r  =  15  a^).  Results  are  calculated  from  (12) 

assuming  f  =  0.15  and  a,  /a  =  0.8.  Since  I  is  positive  but 
°  n  bo  p  n  v 

shows  a  pronounced  dip,  this  beam  exhibits  core  and  halo  behavior 


Fig.  5  Enclosed  currents  1^'  '(r),  I  (r),  and  1^'  '(r)  for  a  case  similar 
to  Fig.  4  but  with  a^/a^  =  0.4.  Since  I^^fr)  reverses  sign  at 
r/a  =  2,  this  case  exhibits  off-axis  hollowing. 


6  Equilibrium  beam  current  density  Jj (r)  as  obtained  by  direct 
integration  of  Ampere's  law  and  the  Boltzmann  equilibrium 
condition.  Also  shown  are  Jp(r),  which  (by  fiat)  has  a  Bennett 
profile  with  radius  ap  and  JBen(,:)>  a  Bennett  profile  fitted  to 
Jb(r)  in  the  vicinity  of  r  =  0.  For  this  case,  T/Tg  =  0.6, 

1-f  =  0.47,  and  there  is  a  weak  halo,  but  overall  JBen(l)  is  a  good 
approximation  to  J,(r). 


1  Cm  rent  densities  J.(i),  J  (i),  and  J  (i)  foi  a  case-  with 

h  p  Ben 


T/Tjj  0.2  and  1  I  0.18.  .J^O  )  is  obtained  in  the'  same  niannn 


as  in  Fig.  (>  and  exhibits  a  pi  ononneed  roie/halo  piolile. 


FiR.  8 


Current  densities  J.(r),  J  (r),  and  Jn  (1)  for  a 
b  p  Ren 

T/Td  0.1  and  1-f  0.083.  This  case  is  siinilai 

D 


case  with 
to  Fig.  7 


the  larger  return  current  results  in  oft  axis  hollowing. 


Fig.  10  Current  densities  J^(r)  (solid  line)  and  J^r)  -  J  (r)  +  J^(r) 

(dashed  line)  taken  from  a  SARLAC  simulation  of  Case  1  at  C  =  SO 

and  7.  -  240  cm.  There  is  no  evidence  for  a  halo  or  hollowing  in 

this  well-matched  (f  =  0.5)  10  kA  case.  The  data  is  taken  alone 

m 

the  positive  x-axis  (9  0).  and  current  densities  are  in 

') 

s ta tamps/cnT  . 


excellent  fit  for  this  case. 


•  t 


13  Plot  of  the  position  and  of  a  subset  of  simulation  particles 
versus  C,  for  the  10  kA  simulation  (Case  1)  shown  in  Figs.  10-12. 

At  7.  =  240  cm,  the  beam  does  not  show  evidence  for  hollowing,  halo 
formation,  or  large  amplitude  hose  motion. 
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15x10 


IOxIO11 


5x  1011 


-5x  1011 


r  (cm) 

Fig.  15  Plot  of  Jb(r)  (solid  curve)  and  Jn(r)  (dashed  curve)  from  a  SARLAC 

simulation  with  I  =  100  kA  and  f  =  0.4  (Case  4).  At  z  =  720  cm 
o  m 

and  C  =  240  cm,  the  beam  forms  a  halo  extending  from  5  to  15  cm 
from  the  axis. 


2 


20 


25 


30 


-5x10 


-9 


10 


15 


r  (cm) 

Fig.  16  Plot  of  pinch  force  F^fr)  vs.  r  for  the  SARLAC  simulation  (Case  A) 
shown  in  Fig.  15.  The  halo  forms  in  the  region  near  the  off-axis 
local  minimum  in  |F  (r)|. 
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y  (cm)  x  (cm) 


100 

50 

0 

-50 

-100 

0  100  200  300  400  500 

f  (cm) 

100 

50 

0 

-50 

-100 

0  100  200  300  400  500 

r  (cm) 

Fig.  17  Plot  of  simulation  particle  positions  x.  and  from  the  core  and 
halo  example  (Case  4)  discussed  above,  taken  at  z  =  720  cm.  The 
dense  core  is  clearly  seen,  and  the  sharp  edge  of  the  halo  at  10-15 
cm  from  the  axis  is  also  apparent. 
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1 


Fig.  18  An  example  of  core  and  halo  behavior  in  a  poorly  matched  (fm  = 

10  kA  beam.  The  plot  of  J^r)  (solid  curve)  and  the  Bennett  fit 
JBen(r)  (dashed  curve)  are  from  Case  2  at  z  =  240  cm  and 
C  =  80  cm.  Most  of  the  current  flows  in  the  halo. 


Fig.  19  An  example  of  off-axis  hollowing  produced  by  a  SARLAC  simulation  of 
a  500  kA  beam  (Case  8).  The  plots  of  J^(r)  (solid  cur\ -)  and  the 
Bennett  fit  Jgen(r)  (dashed  curve)  are  taken  along  the  spoke  at 
0  =  7it/8  which  was  the  direction  of  the  local  maximum  in  J^(r,0). 
The  plot  is  for  C  =  480  cm  and  z  =  1200  cm. 


Fig.  20  Particle  position  plots  at  z  =  1200  of  the  radius  r.  =  (x.  + 

2  1/2 

yj  )  (top  of  frame)  and  y  position  (bottom  frame)  for  the 
off-axis  hollowed  SARLAC  simulation  (Case  8)  described  in  Fig.  19 
A  gap  in  particles  is  apparent  in  the  r.  plot  at  the  top  edge  of 
the  denser  band  at  20-25  cm  from  the  axis.  Geometrical  factors 
obscure  the  gap  when  viewed  from  the  side  (bottom  frame). 


Fig.  23  Beam  displacement  X(z)  taken  from  simulations  of  a  well-matched 
10  kA  beam  (Case  1).  The  three  SARLAC  curves  are  for  C  =  40,  80, 
and  160  cm,  and  the  VIPER  curve  is  C  =  160  cm.  As  expected  for 
this  Bennett-like  example,  the  SARLAC  hose  instability  displacement 
at  (  =  160  cm  agrees  well  with  the  VIPER  result. 


Beam  displacement  X(z)  at  £  =  120,  240,  and  360  cm  taken  from  a 
SARLAC  simulation  of  a  100  kA  beam  (Case  4)  which  exhibits  core  and 
halo  behavior.  Also  shown  is  the  result  from  a  VIPER  simulation  at 
C  =  360  cm.  In  this  comparison,  the  SARLAC  model  predicts  much 
more  hose  growth  than  does  VIPER. 


cm 


100 


x  -  f  cross-section 


£  (cm) 

Fig.  26  Simulation  particle  positions  for  the  core  and  halo  example  (Case 
4)  taken  at  z  =  1200  cm.  The  hose  instability  distortions  of  the 
beam  core  are  clearly  observable.  (In  Fig.  17,  which  was  taken 
from  the  same  simulation  at  a  shorter  propagation  distance,  the 
instability  had  not  yet  grown  into  the  large  amplitude  regime). 
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Fig.  27  Beam  radius  ^or  the  core  an(^  halo  simulation  (Case  4)  shown 

in  Figs.  25  and  26.  At  z  =  720  cm,  the  half  current  radius  from 
the  SARLAC  code  quickly  reaches  a  nearly  constant  value  of 
10-12  cm,  although  the  core  radius  is  much  smaller  (~  aQ  =  3  cm). 
The  VIPER  code  at  z  =  720  cm  predicts  a  beam  radius  which  tapers 
significantly. 
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